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A method for selecting optimum NASTRAN analysis set degrees of freedom for the 
dynamic eigenvalue problem is described. Theoretical development of the Guyan 
reduction procedure on which the method is based is first summarized. The 
algorithm used to select the analysis set degrees of freedom is then developed. Two 
example problems are provided to demonstrate the accuracy of the algorithm. 


10 INTRODUCTION 

A NASTRAN user is faced with two major difficulties when solving a dynamic 
eigenvalue problem. First, an eigenvalue solution is expensive for most structural 
problems encountered in engineering applications. Second, many more degrees ot 
freedom (DOF) are required to define a structure's elastic properties than are 
required to define its inertial properties, which tends to exacerbate the first 
difficulty. 


A popular method for easing the severity of these difficulties is to reduce the 
problem size using Guyan reduction (Reference 1). This method allows the user to 
preserve the elastic properties of the reduced problem set while reducing the 
problem size to one more manageable for a dynamic eigenvalue solution. At the same 
time, the mass properties are also condensed with some penalty associated with the 
reduction of mass from the coordinates being eliminated. The present paper 
describes an approach for optimizing the partitioning process to minimize this 

penalty. 

Theoretical development of the Guyan reduction method is presented in Section 2. 
Section 3 describes the algorithm used to select automatically the analysis set degrees 
of freedom. Verification of the method is presented in Section 4. Conclusions arc 
presented in Section 5. 


THE GUYAN REDUCTION METHOD 


By way of introduction, the Guyan reduction method will first be reviewed. 


The dynamic eigenvalue problem is given by the equation 

([K] - MM ]){«>} =0 O) 


where K is the structural stiffness matrix, M is the structural mass matrix, X is the 
eigenvalue, and <J> is the eigenvector or modal displacements. The Guyan reduction 
method starts by partitioning Equation 1 into independent DOF, designated in 
NASTRAN as the A-set, and dependent DOF, designated as the 0-(for OMIT) set. After 
performing this operation Equation 1 becomes 


175 


( 2 ) 


/r 


Kaa K 


ao 


VL K ao K oo J 


^aa ^®ao 
L M ao M ooJ 


to- 


J 


where the subscript "a" denotes A-set DOF and the subscript "o" denotes O-set DOF. 

A set of constraints for the O-set displacements can be derived by solving for $ 0 in 
terms of 4> a using statics, or 

K ao*a + K oo*o = 0 

The O-set displacements now become 

♦o = G o^a (4) 


where 


Go = -l<iX> (5 > 

Equation 4 defines 4» 0 as the deflections at O-set DOF due to unit displacements at the 
A-set DOF. Stated another way, the O-set displacements, 4> 0 . arc constrained to move in 
relation to A-set displacements, $ a , as governed by the transformation matrix G 0 . This 
relationship constitutes a Ritz transformation of the eigenvalue problem. The 
transformation written in terms of the full displacement set is 
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Using this Ritz transformation the reduced mass and stiffness matrices become 

|MJ = [G] T [M][G] (7) 


and 

[KJ = [G] T [K][G] 


( 8 ) 


Performing these operations on the matrices in Equation 2 we get 


|M m I = IMJ + 2[M a0 ][G 0 ] + [G„] T [M (>0 )[G () ] 


(9) 


and 


= [K aa ] + [K a0 ][G 0 ] 


(10) 
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The mass of the system will be redistributed based on the elastic connections between 
the O-set DOF and the A-set DOF as shown in Equation 9. 


Note that Guyan reduction is exact when Moo (®nd hence M a o) ts a null matrix and 
gives the best solution for any selected partition when it is not. It does not, however, 
address directly the problem of selecting most effectively the set of independent DOF 
that will best serve the aims of the user. For this, a means of removing terms from 
the mass matrix so as to minimize the impact on the solution accuracy must be 
determined. 

3 0 ASET SELECTION ALGORITHM 

As stated previously, Guyan reduction is exact when M 0 o * s nu H< or when the O-set 
mass to stiffness "ratio" is zero. As the mass to stiffness "ratio" between Mq 0 and Kqo 
increases, the accuracy of the Guyan reduction method decreases. This 
generalization forms the basis of the A-sct selection method. 

The six step method for determining the A-set DOF is as follows: 

1. Execute NASTRAN to obtain an initial M aa . K aa , and A-set table. The mass and 
stiffness matrices can be reduced as desired in NASTRAN as long as the modal 
content over the frequency range of interest is retained. Note that no reduction 
need be performed at this stage but the initial constraint equation must be 
applied. 

2. Define the number of DOF that will be in the final A-set. These DOF may also 
contain a "kernel" set of DOF that will remain in the A-sel regardless of their 
mass to stiffness ratio. 

3. Determine the minimum mass to stiffness "ratio" for the O-set DOF. Because M 
and K are diagonally dominant, this ratio is most easily approximated by 
stripping the diagonal from M and K and scanning for the minimum Mij/Kjj 
which we will call min(M/K). The min(M/K) DOF is then partitioned from M and 
K and reduced from the system, provided it is not a member of the kernel set. 

4. Repeat step 3 until the desired number of DOF remain in the A-set. 

5. Write NASTRAN ASET bulk data cards for the retained DOF 

6. Check the A-set to determine if desired modes are adequately defined. 


To improve the efficiency of the check process, the mass and stiffness matrices may 
be saved during Step 5. These matrices can then be used in an eigenvalue analysis to 
determine if the selected A-sct is adequate. 

The user may, if desired, decide to refine the A-set further if it is concluded that more 
DOF can be reduced from the problem. To simplify this second reduction, the A-set 
listing and matrices from Step 5 can be used as input to Step 2. The process would 
then proceed as before. 

Occasionally, too few DOF will be defined in the A-set. By keeping track of the DOF 
placed in the O-set during each iteration, the user may simply review DOF that were 
omitted during previous iterations to determine DOF that are required to define the 
mode or modes lost because of the Guyan reduction. He may then selectively include 
those DOF deemed necessary to the A-set by adding these DOF on his ASET bulk data 
cards. Alternatively, he may save intermediate ASET card images for convenience. 
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Because the algorithm currently works on one DOF at a time, the user should use 
NASTRAN to make the problem size as small as possible to decrease the solution time. 
Though reducing several DOF during each iteration is a desirable feature, no 

definitive method for including this feature in the algorithm has yet been developed. 

More information on this topic is presented in the conclusions. 

The algorithm described above virtually guarantees that the smallest A-set will be 
obtained with minimal effort, provided that too severe a reduction is not specified. 

The general procedure for selecting the A-set automatically should be clear from the 
discussion above. The process is best illustrated, however, by performing sample 

calculations on a simplified model, as shown in the next section. 

4 0 METHOD VERIFICATION 

Two sample problems were developed to validate the A-set selection method. The first 
problem is a simplified model of a three story building. The reduction operations are 
performed by hand to clarify the algorithm. The second sample problem determines 
the A-set of a 3600 DOF NASTRAN model. The A-set for this problem was generated 
using a program developed by McDonnell Douglas Space Systems Company-Huntsville 

Division (MDSSC-HSV). The data from these sample problems verify the algorithm 
outlined in Section 3. 

The simplified model of the three story structure is shown in Figure 1. The mass and 
stiffness matrices are also shown. The fundamental frequency of this system is 
1.45 Hz. We want to reduce the problem to a one DOF system. 
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Figure 1. Simplified Three Story Building 
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First wc find the min(M/K) for this system which is 2/2000 = 0.001 for displacement 
U 3 . Partitioning this DOF from M and K yields 
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The G 0 matrix for this problem is 

g .= 255 o «° - 80 °1 = [00 oa] 

The reduced mass and stiffness matrices are found using Equations 9 and 10 and are 
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We repeat the steps to determine the mass and stiffness of the one DOF system. 
Performing these steps produces M = 2.48 and K = 218.2. The frequency for this one 
DOF system is f\ - 1.50 Hz which is 3.5 percent higher than the exact frequency of 

1.45 Hz. 


Though the frequencies show excellent agreement, correlation between the mode 
shapes should also be verified. Back transforming using G 0 we get 


{<M = 
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for the one DOF system. We will use the modal assurance criterion (MAC) described in 
Reference 2 to measure the correlation between this mode shape and the exact 
mode shape. The MAC between any two modes varies from zero, meaning no 
correlation, to unity, meaning perfect correlation. The MAC for these modes is 0.987 
indicating that little modal accuracy was lost during the reduction. 

The second sample problem involves finding an A-set for the model shown in 
Figure 2. The unreduced model has approximately 3600 DOF. Currently, the model 
A-set has 180 DOF which was used as a starting point for this problem. This A-set was 
funher reduced to 50 DOF using the MDSSC-HSV developed program based on the 
selection algorithm described in Section 3. The final A-set size is approximately 25 
percent of the original A-set size. 
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Figure 2. NASTRAN Model for Sample Problem 2 


Table 1 shows a comparison between the frequencies and mode shapes of the 180 DOF 
model and the SO DOF model. The frequencies show excellent agreement with a 
maximum difference of 1.4 percent for the sixth mode. The mode shapes are almost 
perfectly correlated between the the 180 DOF model and the 50 DOF model. Indeed, it 
may be possible to reduce the problem size even further. 


Table 1. Frequency and Mode Shape Comparison 
Between 180 DOF Model and 50 DOF model 

MODE 

/l 80 

/50 


MAC 

1 

11.9 

11.9 


>0.999 

2 

12.9 

13.0 

0.8 

>0.999 

3 

24.1 

24.2 

0.4 

0.998 

4 

24.9 

25.0 

0.4 

0.996 

5 

33.1 

33.3 

0.6 

0.998 I 

6 

62.3 

63.2 

1.4 

0.992 


5.0 CONCLUSIONS 

A method for automatically selecting the NASTRAN A-set DOF was described. 
Theoretical development and an outline of the steps involved were provided. Two 
example problems were provided that demonstrate the use and the accuracy of the 
method. Some potential enhancements have been identified and will be briefly 
summarized here. 



































One ootential enhancement noted earlier would be to reduce multiple DOF during 
each iteration Because of the redistribution of the mass of the system, simply 
Sffn,* percentage of Ute DOF a. each itemion is to be discouraged. The 

reason for this is best demonstrated with an example. 

Consider the simply supported beam of Fipire 3. Because alii rf the DOF have identical 

mass to stiffness ratios, the removal would begin with the first DOF 

min(M/K) If a 20 percent reduction rate were chosen then uj and U2 

remold in the first iteration, which could ultimately result in a poorly chosen A-set. 
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Figure 3. Simply Supported Beam 


A second potential enhancement would be including a method in the > 1 ** 

would determine the optimum number of A-set DOF based on a user defined uppe 
bound frequency of interest. Because the algorithm removes terms with a high 
pseudo frequency, i.e. large Kji/Mn, an approach based on the pseudo frequencies of 
the reduced system could be used to predict the minimum required number of A-set 
DOF. 

Even without these enhancements, the method has been successfully implemented at 
MDSSC-HSV. The often tedious, and sometimes eiror prone A-set selection process nas 
been automated, saving engineering time while increasing A-set efficiency. 
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